Bose-Einstein condensation of 2D dipolar excitons: Quantum Monte Carlo simulation 
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The Bose condensation of 2D dipolar excitons in quantum wells is numerically studied by the dif- 
fusion Monte Carlo simulation method. The correlation, microscopic, thermodynamic, and spectral 
characteristics are calculated. It is shown that, in structures of coupled quantum wells, in which 
low-temperature features of exciton luminescence have presently been observed, dipolar excitons 
form a strongly correlated system. Their Bose condensation can experimentally be achieved much 
easily than for ideal or weakly correlated excitons. 
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I. INTRODUCTION 

Two-dimensional (2D) dipolar excitons, in which elec- 
trons (e) and holes (h) are spatially separated, can be 
created in structures of coupled (CQWs) quantum wells 
(QWs) and single QW (SQW) in a strong external 
transverse electric field Q. At low temperatures, e and 
h in these systems reside in the lowest transverse (with 
respect to the QW plane) quantized states. In this case, 
the dipole moments of excitons (in the ground axially 
symmetric state) are all directed perpendicular to this 
plane. There are two types of excitons in CQWs: direct 
excitons, whose electrons and holes are in the same QW, 
and spatially indirect excitons, whose electrons and holes 
are located in different QWs (see details in experimental 
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second case, exciton dipole moment is d — eD, where 
the quantity D is determined by the geometry and is 
equal to the characteristic interwell distance (see Eq. © 
below). In single quantum wells, a dipolar exciton can 
occur due to a strong transverse electric field, which in- 
duces a dipole mome nt B , EH- In both cases, unlike 
quasi-2D atomic gases [l6| , two 2D dipolar excitons can- 
not pass "over one another" in the QW plane. Dipolar 
interaction of 2D excitons is theoretically studied in Rcf. 

At low exciton densities, 2D dipolar excitons are far 
from each other, so exchange effects between electrons 
(holes) of neighboring excitons are suppressed by a dipole 
barrier, related to a exciton dipole-dipole repulsion, 



P 



exp 



y/m(V(r) - 2E /N)dr 



Here, P is the probability for a 2D dipolar exciton to 
tunnel through the dipole barrier; m is the exciton mass; 
Eq is the ground-state energy of the excitons (their ki- 
netic and dipolar-interaction energies); N is the exci- 
ton number; V(r) is the dipolar exciton potential; a ex 



is the exciton diameter (which is of the order of the av- 
erage distance a e h between e and h in the exciton); and 
a c i is the classical turning point for the potential V(r) 
(V(a c i) — 2Eo/N). In this case, 2D dipolar excitons can 
be considered as (composite) bosons. At high exciton 
densities, at which the distance between neighboring ex- 
citons is smaller or of the order of the exciton radius, the 
Fermi properties of the internal exciton structure play 
a significant role. This leads to a transformation of a 
dipolar exciton gas into an e-h system lj| ]j| (although, 
in the isotropic system at T = 0, the crossover from an 
exciton Bose condensate to a BCS state, with the gap 
decreasing with an increase in the density, occurs [H). 

The possibility of the Bose condensation of 3D exci- 
tons was under consideration long ago [2p| . However, in 
a homogeneous infinite 2D system, because of the ultra- 
violet divergence of the exciton momentum distribution, 



the Bose condensation is possible only at T = [2l| . In 
the range of finite temperatures, only a quasi- long-range 
(power-law) order [22I ] and a superfluid phase can occur, 
which are destroyed due to the dissociation of pairs of op- 
positely directed vortices and formation of free vortices 
at the Berezinskii-Kosterlitz-Thouless (BKT) transition 



point [23J, [2J, |25|, |26|. However, in finite 2D systems 



(traps), the ultraviolet divergence of the momentum dis- 
tribution is cut off, so, the Bose condensate can exist 
fliL [2^ . This opens opportunities for studying the Bose 
condensation of 2D dipolar excitons in real experimental 
systems. 

At present, the search for the Bose condensation of 
2D dipolar excitons in CQWs 0, 0, !, i] and in SQW 
has attracted much attention from researchers. Inter- 
esting optical effects observed in studying the lumines- 
cence of 2D dipolar excitons in CQWs at low temper- 
atures point to significant experimental progress in the 
investigation of the collective properties of 2D dipolar ex- 



citons @, S, 0, El EI EI EE BJ. The observations of 
a narrow-beam luminescence from a 2D dipolar exciton 
Bose condensate along a normal to the QW plane and 
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of the off-diagonal long-range order of Bose condensed 
2D dipolar excitons in the coherence experiment have re- 
cently been performed by Timofeev et al. in SQW in a 



transverse electric field 15 



Theoretically, for the collective state of 2D dipolar ex- 
citons in QWs, superfluidity 28, 2(J[3(|, Bose condensa- 
tio n [181 l30l . l3l| (in the exciton regime) , a BCS-like state 
0, l32j (in the e-h s yste m) strong correlations @, (33|, 
crystallization 0, Hi , H, 3y, 37 1 and a mesoscopic super- 
solid pha se [iH , kinetic and relaxation processes , vor- 
tical |40( and quasi-Josephson 41| effects, optical proper- 
ties [42j as well as these in circular traps [43j , a superfluid 
transition (crossover) in traps [3], as well as the effect 
of random fields 2iJ 45|, were studied. 

The superfluidity of Bose condensed dipolar excitons 
was also studied in two-layer e-e systems in a strong mag- 
netic field with the half-occupancy of the Landau levels in 
each layer (v = 1/2 + 1/2; see, e.g., the theory in Ref. [3| 
and experiments in Ref. [47j and references therein) . 

2D dipolar bosons were also studied in rarefied atomic 
gases 48 1 and polar molecules (iij . 

Among the above studies into the exciton Bose con- 
densate, the model of an ideal exciton gas and that of a 
weakly correlated exciton gas were predominantly used as 
theoretical approaches. Unfortunately, these approaches 
have a rather limited domain of applicability. Indeed, 
at low densities, the repulsive dipole potential can be 
described by only one parameter, the length a s of the 
isotropic s-scattering. Therefore, the properties of such 
a potential should be universal, i.e., the same for all in- 
teraction potentials having the same scattering length; in 
particular, the properties should be the same as for the 
system of hard disks with the diameter a s . However, the 
latter system is known to be weakly correlated only in 
the ultrararefied case [H(J. Consequently, the model of 
weakly correlated excitons is valid only for ultrararefied 
systems, whose critical temperature is extremely low. For 
real experimentally observed densities, this model can 
give only a qualitative description. Therefore, the quan- 
titative theoretical investigation of the 2D QW dipolar 
exciton collective state requires a more accurate model 
and more precise, desirably, ab initio, calculations. 

This work is devoted to the detailed microscopic in- 
vestigation of 2D dipolar excitons by means of the nu- 
merical simulation by the diffusion Monte Carlo method. 
The main result of our simulation is the inference that, 
in the experimentally observed low-temperature regime 
d 11 0, 2D dipolar excitons in CQWs form a 
strongly correlated system [Hl|. We find the following 
evidences for this fact. 

(i) The dimcnsionless adiabatic compressibility £ = 
m 3 / (2ith 2 x) (x is the corresponding dimensional quan- 
tity) and the contribution f = /iim/(27rfi, 2 n) to the 
chemical potential fii related to the dipolar exciton- 
exciton interaction prove to be much greater than unity, 
C, (' ^> 1 (Here, n is the total exciton density in all the 



g ex spin degrees; g ex — 4 for GaAs). For weakly corre- 
lated excitons, £, (' <C 1. 

(ii) The density no of the exciton Bose condensate (in 
all the g ex spin degrees) at T = proves to be two to four 
times smaller than the total density n (in the regime of 
weak correlations, uq w n). 

(iii) A clearly pronounced hump is observed both in 
the pair distribution function and in the structure factor, 
which points to the presence of the short-range order. In 
addition, at higher exciton densities, e.g., at n = 5 • 10 10 
cm~ 2 , two humps are observed in the structure factor 
and three in the pair distribution function. 

(iv) The spectrum of excitons is very far from the 
Bogolyubov shape (although, even at the boundary be- 
tween the exciton and e-h regimes, the roton minimum 
is not yet reached under the experimental conditions of 
Refs. 0,0, BBS). 

(v) The local superfluid exciton density at the super- 
fluid transition temperature, ni(T c ), is close to the total 
exciton density n (whereas, for weakly correlated sys- 
tems, the quan tity ni (T r ) is logarithmically small com- 
pared to n [53J) [541 ]. 

(vi) The temperature of the superfluid transition in 
the corresponding infinite system is only slightly less 
than the degeneracy temperature T^eg, while the quasi- 
condensation temperature is 2-2.5 times higher than Td eg 
(for the spin-depolarized exciton gas in all the g ex = 4 
spin degrees). In contrast, according to the theory of 
weakly correlated excitons, these temperatures are loga- 
rithmically small compared to Td eg j53| . 

(vii) The profile of the exciton Bose condensate in a 
large 2D harmonic trap at T — appreciably differs from 
the Thomas-Fermi inverted parabola. 

In some structures of CQWs ToL ITl [ and SQW 15 1, the 
condensed gas of 2D dipolar excitons forms a system with 
intermediate correlations (£, £' ~ 1). In this case, the 
above evidences weaken (evidence (iii) vanishes). How- 
ever, in structures of CQWs and SQW with the width 
of these wells being sufficiently large, at high densities of 
excitons, their correlations prove to be so strong that a 
roton minimum appears (and even the crystallization of 



36|,|37|). 



excitons is possible; see, e.g., Refs. 

The paper is organized as follows. In Section [Til the 
model is discussed, the problem is formulated, and the 
relation with experiment is considered. Also, the details 
of simulation are presented. Section IIIII is devoted to a 
homogeneous exciton system. A detailed analytical pro- 
cessing is constructed based on the ab initio numerical 
simulation, with the microscopic parameters of the prob- 
lem being calculated. Using these results, in Section \W\ 
we analytically investigate a large harmonic trap in the 
local density approximation. The possibilities of exper- 
imental observation of strong correlations in the system 
of 2D dipolar excitons in CQWs are discussed in Section 



Vj In section [VI] we conclude the study. 
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II. MODEL AND CALCULATION METHOD 

Two 2D spatially separated excitons repel each other 
as dipoles, V(r) oc 1/r 3 , if the distance between them is 
much longer than the effective spacing D between e and h 
layers. In the exciton regime, the distance between neigh- 
boring excitons in structures of Refs. [1, @, 0, H B 0, U 
ll2l . ll3L 0, EH , according to a rough estimate (see details 
in [55(), exceeds D. Therefore, the model of dipolar exci- 
tons is a good approximation for experimentally studied 
systems. The suppression of the exchange interaction by 
the dipole barrier in the exciton regime makes it possible 
to consider 2D dipolar excitons as structureless bosons 

As a result, at r 3> D, we arrive at the following Hamil- 
tonian of the homogeneous system of structureless 2D 
dipolar excitons: 



Here. 



H = 



x 



O + r 3 . 
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(1) 



(2) 



AitEe^h 2 eh 2 

is the parameter with the dimension of length for the 
exciton dipole-dipole potential V(r), 



V(r) 



h 2 xq 



Xq 



Xq 



8 2 /dx 2 + d 2 /dy 2 + d 2 /dz 2 



d = eD is the dipole moment of an exciton, e > is the 
hole charge, e is the dielectric constant of the structure, 
Sq = 1 /4-7T is the dielectric permittivity of vacuum, and 



D 



(\Mz)\ 2 ~ \Mz)\ 2 )zdz 



(3) 



is the effective distance between the e and h layers, where 
ipe(z) and "4>h{z) are, respectively, the electron and the 
hole wave functions along the z axis. 

We impose periodic boundary conditions. This, on the 
one hand, makes the system finite and, on the other hand, 
imitates its homogeneity. In this case, the momentum 
becomes discrete but remains good quantum number. 

At the first stage, in calculations by the variational 
Monte Carlo method, we specify the trial function of the 
ground state of excitons to be a Bijl-Jastrow function 



(4) 



i<j 



where, 



AKb(2/Vr), r<f c , 
Bexp(-C/f-C/(L-f)), f c <f<L/2. 




(5) 



2 3 
rVn 

Figure 1: Trial wave function f(r) at different exciton den- 
sities n. Curve i (i = 1,2. ..12) corresponds to the density 
nx\ = 2 1-9 . Circles show the point r c of joining of two asymp- 
totics. 



Here, K q(x) is the modified Bessel function of the sec- 
ond kind (the Macdonald function) of the zeroth order, 
r = r/xo, and L — L/xo is the dimensionless size of the 
system L. The constants A, B, and C are chosen from 
the condition f(L/2) = 1 and continuity conditions of 
the zeroth and first derivatives at f = f c , and the joining 
point f c is used as the (only) variational parameter. We 
expect that the trial wave function (U) is a good approx- 
imation for the wave function of the ground state. When 
two excitons approach close each other, the influence of 
remaining exciton pairs becomes small, and the function 
/(f) corrresponds to the solution of the two-particle scat- 
tering problem from the dipole-dipole potential. At short 
distances, function ([5]) is chosen such that it corresponds 
to the exact solution of the scattering problem for zero 
energy. We also verified that the use of the numerical so- 
lution to the scattering problem at finite energy as a two- 
particle factor of /(f) in the Bijl-Jastrow function weakly 
varies the variational energy of the system and almost 
does not affect the diffusion value of the energy. The 
long- wavelength behavior of the function /(f) is deter- 
mined by the collective properties and corresponds to the 
phonon propagation in the 2D system: /(f) cx e - const / r 
[561 ] . The second term is added to the exponent of the 
exponential in Eq. ([5]) in order to satisfy the condition 
that the derivative at the boundary of the box is zero: 
f(L/2) = 0. The function /(f) is shown in Fig. 1. 

As is seen from Eqs. (JlJ, (j4j, and ([5]) the Hamiltonian 
and the trial function are reduced to the dimensionless 
forms that are expressed, respectively, in terms of the 
linear units xq and energy units jmx\. As a result, 
the dimensionless exciton density n = nx\ is the only 
dimensionless parameter of the problem. 

The quantity 1/xq yields the 2D exciton density n, 
which corresponds to the dimensionless density n = 
tiXq = 1 and which can be calculated by Eq. @ . Now, we 
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will consider the range of variation of the dimensionless 
density n in various real experimental systems. 

(1) For the structure of GaAs/AlAs CQWs studied 
by Timofeev and colleagues HQ (D = 13.6 nm, m = 
0.22m (mo = 9.1083 • 10~ 28 g is the free electron mass), 
and e = 12.5), \jx\ = 2.64 • 10 10 cm- 2 . 

(2) For the structure of GaAs/AlGaAs CQWs (Butov 
et al. dHI, D = 12.3 nm, m = 0.22m , and e = 12.5), 
1/xl = 3.94- 10 10 cm" 2 . 

(3) For the structure of asymmetric GaAs/AlGaAs 
CQWs (Moskalenko et al. [§],£> = 14.1 nm, m = 
0.22m , and e = 12.5), 1/xl = 2 - 28 ' 1Ql0 cm ~ 2 - 

(4) For the structure of InGaAs/GaAs CQWs (Snoke 
et al. 03, El], D = 10 - 5 nm > m = 0.14m , and e = 12.5), 
1/xg = 1.83- 10 11 cm- 2 . 

(5) For the structure of GaAs/AlGaAs CQWs (Snoke 

et al. [2E3, D = 15 ' 5 nm ' m = °- 22m o, and £ = 12 -5), 
1/xl = 1.56- 10 10 cm" 2 . 

(6) For the structure of GaAs SQW (Timofeev et al. 

[l5j, Dw6nm,m= 0.22m , and e = 12.5), n = l/xg ~ 
7 • 10 11 cm -2 , which corresponds to the dimensionless 
density n = 1. 

Next, we have to know the range of variation of the 
2D exciton density n in their experimentally achievable 
condensed state. 

The upper limit of the density is the boundary between 
the strong and weak coupling regimes (corresponding 
to the condensed exciton state and to the BCS-like e-h 
state). For structure (1) [B|, Hi with an average electron- 
hole distance of a e h = 17 nm [7|), the corresponding crit- 
ical de nsity determined according to the rough estimate 
of Ref. 



55J is 



n max = 0,11 2 ge:r ~ 5.17-10 10 cm 2 , n max « 1.96 (6) 

(here, it is taken into account that, in GaAs, excitons are 
in the g ex — 4 spin degrees). 

Therefore, characteristic exciton densities that can be 
achieved in modern experiments in QWs, at which low- 
temperature features in exciton luminescence are ob- 
served d @, 0, 1, S], lie in the range 10 10 < n < 5 • 10 10 
cm- 2 , i.e., 1/4 < n < 4. 

However, in very wide (wider than 100-200 nm) 
SQWs based on GaAs, such high dimensionless densi- 
ties (n max > 290 j3?| ) can be achieved at which excitons 
crystallize [36j ]. 

We performed simulation in the range of the dimen- 
sionless density 1/256 < n < 8. For structure (5) [l2l[l3j 
(in which the exciton correlations are stronger than in 
other structures), at the density n = 8, the dimensional 
exciton density is high: n = 1.2 ■ 10 11 cm -2 . In struc- 
ture (6) [HI (where the exciton correlations are weaker 
than in other structures), at the dimensionless density 
fi = 1/256 (and at the corresponding dimensional density 
n w 2.7- 10 9 cm -2 ), the superfluid crossover temperature 



of excitons in a trap of about 100 /mi in size is low: 0.12 
K (see Q and Section [nTF| . 

We should note that, in a sufficiently weak electric 
field, very small dipole moments eD and very low di- 



mensionless densities n 



HXq 



cx D 4 



Eq. 0) 



be realized in SQWs. However, we are not interested 
in this case. Indeed, at D — > 0, the contribution of 
the van der Waals potential to the exciton scattering 



length a s (a 



WDW 



a e h) considerably exceeds the con- 



tribution from the dipole potential to this parameter 
(af = 3.17222xo cx D 2 — s- 0). As a result, excitons 
efficiently loose the property of dipolarity. 

Nevertheless, in GaAs-based SQWs (m = 0.22m , e = 
12.5, a e h « 12 nm), in a sufficiently weak electric field, 
the dimensionless density of a rarefied 2D dipolar exciton 
gas with a superfluid crossover temperature of about 0.1 
K and with a dipole scattering length of the order of a e h 
can be rather low: n - 2" 11 = 1/2048. 

We performed the simulations by the quantum Monte 
Carlo (MC) method for 12 different densities n = 2 l ~ 9 
(i = 1, 2. ..12) in the range 1/256 < n < 8. The number 
of excitons was taken to be TV = 100. Initially, we per- 
formed the calculations by the variational MC (VMC) 
method with the trial function defined by ([4j) , ([5j) and 
optimized the parameter f c . At the second stage, the 
calculations were performed using the ab initio diffusion 
MC (DMC) method [57] . in which, to accelerate the con- 
vergence, the same trial function was taken but with the 
parameter f c being already optimized. 

As is known, the inaccuracy of the trial function intro- 
duces an additional error in the calculation of correlators 
by the DMC method. To reduce this error, the data 
obtained by the two variants of the MC method were 
linearly extrapolated, which allowed us to eliminate the 
first order with respect to the small difference between 
the exact and trial functions of the ground state [H^ . 



III. MEASUREMENT RESULTS. A 
HOMOGENEOUS SYSTEM 

In sections IIIII and IIV1 we will use the dimensionless 
system of units h = m = Xq = 1, so, n = n, L = L, etc. 



A. The the ground state energy. The chemical 
potential. The adiabatic compressibility 

The total energy of the ground state of a two- 
dimensional dipolar exciton can be represented as the 
sum of two terms, 



E,, 



2irNne 2 D 



E . 



(7) 



The first term corresponds to the energy of a capacitor 
formed by plane-parallel e and h layers spaced by the 
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0.01 0.1 1 10 

density 

Figure 2: (a) Ground state energy per one exciton, Eo/N, 
as a function of the density (squares are the DMC data) and 
polynomial fitting curve (|lip (we set h = m = xo = 1). (b) 
Dimensionless characteristics analytically calculated from the 
polynomial fitting of (JTTJ as functions of the density: energy 
C per one exciton; contribution £' to the chemical potential 
related to dipolar interactions; and adiabatic compressibility 
C at T = 0. 



effective distance £). The second term is responsible for 
the kinetic energy of N dipoles and for their dipolar in- 
teractions. It is equal to the energy of the ground state 
of the system of these dipoles. This energy is calculated 
by the DMC method. 

The exciton chemical potential counted from the exci- 
ton band edge can be expressed via their total energy 



dE t , 



dN 



Me + Mi 



(8) 



where fj, e = Ai:ne 2 D/e is the electrostatic contribution, 
and 



Mi 



dEo 
dN 



(9) 



is the contribution to the chemical potential of exci- 
tons, which corresponds to their interactions. This 
term defines the adiabatic compressibility of excitons 
X, 1/x = J(6fJ,i/Sn(r))dr. For a homogeneous system 
(n(r) = n), the variational derivative is reduced to the 
ordinary derivative (J{Sfii/Sn(r))dr = dfii/dn), so 



1/X 



djjj, 
dn 



(10) 



Figure 2(a) presents the results of the DMC simulation 
of the ground-state energy of interacting dipoles per one 



particle, Eo/N, at different densities. The calculated val- 
ues lie with a high accuracy (within the limits of 0.025%!) 
on the curve of the polynomial fitting 

Eq/N = a e exp(6 e \nn + c e In 2 n + d e In 3 n + e e In n), 

(11) 

where a e = 9.218, b e = 1.35999, c e = 0.011225, d e = 
-0.00036, and e e = -0.0000281. 

Figure 2b shows the energy per particle, Eq/N, the 
contribution to the chemical potential [ii that is caused 
by dipolar interaction, and the adiabatic compressibility 
X which were analytically calculated from the fitting to 
energy pT|) (see Eqs. ([9]), (fTT)|) ) and were expressed in 
terms of the quantities and £ respectively, 



E /N = wnC 



1/X = 2tt(. 



(12) 



(13) 



(14) 



We find from this simulation that, in the entire range 
of calculated densities, the dimensionless adiabatic com- 
pressibility (i.e., the dimensionless interaction) proves to 
be greater than or of the order of unity, 



C> 0.774- 1, l/256 <n< 



(15) 



In the regime of high densities (1/4 < n < 4), we 
find that ( 3> 1 (see Fig. 2(b)). This testifies to strong 
correlations in the dense gas of 2D dipolar excitons in 
CQWs, whose condensation was studied in Refs. (3,0,0, 



At low densities, the quantity £ *~ 1. This is the 
regime of intermediate correlations, which was realized 



in Refs. [Id, El El 



The regime of weak correlations ((< 1) has not been 
realized even at the lowest calculated density n = 1/256 
(at which the superfluid crossover temperature in struc- 
tures considered in Refs. i, i, 0, B i, E3, El El El, Ei, 
E3| is not higher than 0.12 K; see Section HI]) . 



B. The one-body density matrix. The Bose 
condensate fraction. The microscopic phonon length 

scale 

On long-wavelength scales r that exceed a certain value 
ro (i.e., at r 3> ro), for the polar-circle- averaged one-body 
density matrix 



/■27T 

pi(r) - / (*+(r)*(0))^/27r 
Jo 



(16) 



of Bose condensed 2D dipolar excitons, the following hy- 
drodynamic expression holds at T = (see Appendix 



6 



A): 
P\\r) 



1.0 



2tt / 

no [explj2e ipI 
o \p^° 



2Np 2 



AN 




r > ro- 

In Eqs. (fl"6]l and (fT7|) ^(r) is the exciton field operator, 
(...) is the averaging over the exciton ground state, tp 
is the polar angle of the vector r, e p is the excitation 
spectrum at T = (which we approximately determine 
from the structure factor by the Feynman formula; see 
Sections IIII D[ IIII Ep , and ro is the characteristic micro- 
scopic phonon length scale of the system that separates 
the long-wavelength (hydrodynamic, r 3> ro) and mi- 
croscopic (r <C ro) ranges (Fig. 3(b)) and corresponds 
to the ultraviolet limit of applicability of the hydrody- 
namic method. The fact that the sum in Eq. (fT7)) is dis- 
crete with respect to momentum (p = (27r/L)n, where 
n = {n x ,n y } ^ is the integer- valued vector) allows us 
to correctly take into account both the finiteness of the 
size of the system and the periodic boundary conditions. 

The MC data for the one-body density matrix pi(r) 
(the linear extrapolation of the VMC and DMC data) at 
different exciton densities are shown in Fig. 3(a). 

As shown in Fig. 3(b) the excelent coinciding of VMC 
and DMC results (within the line thickness) is the ev- 
idence of a good choice of the trial wavefunction and, 
consequently, a high accuracy of our simulation. 

Figure 3(e) shows the one-body density matrix pi(r) 
calculated at n = 1 and large r. This calculation ideally 
coincides (within the line thickness) with the hydrody- 
namic equation (fT7|) . This testifies to the internal con- 
sistency of our MC simulation, to the high measurement 
accuracy, and to the validity of the hydrodynamic de- 
scription in this range. 

For comparison, Fig. 3(e) presents the asymptotic form 
for the one-body density matrix at large r in the corre- 
sponding infinite system (L — > oo) at n = 1, 



1 



Cs 1 

47m r 



L — > oo, r > r , l/c s 



(18) 

which was obtained from Eq. (|17|) by the formal replace- 
ment of the sum over p by the integral. In Eq. (Ti~g|) , the 
exponential is expanded into a series, c s — Wn/x 1S the 
sound velocity at T = 0, and njf = 0.3405 ± 0.001 is the 
density of the Bose condensate in the infinite system at 
n=l. The quantity ng was obtained by quadratic ex- 
trapolation of the values of p\(L/2) to the macroscopic 
limit with respect to the powers of 1/iV using 14 values 
of the total exciton number N from 25 to 200. 

In Fig. 3(c), the density dependence of the the exci- 
ton Bose condensate fraction is presented, which was ob- 
tained by the best fitting of hydrodynamic equation (|17D 
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Figure 3: (a) One-body density matrix pi(r) at T = for 
12 different densities (n = 2 i_9 , i = 1,2. ..12). (b) One-body 
density matrix calculated by the VMC and DMC methods, as 
well as linear extrapolation (LE) of VMC and DMC results 
(thick curves) and hydrodynamic fits (thin curves) at n = 1 
(the curves practically coincide) . The arrow shows the charac- 
teristic microscopic phonon length scale ro (in units of 1/ ^/n) . 
(c) Bose condensate fraction at T = (n/no) and polynomial 
fitting by Eq. (19) as functions of the density, (d) Character- 
istic microscopic length scales of the system: phonon scale ro , 
ideal-gas scale (see Section [III E|) ric, healing length £, and 
interexciton distance 1/\A* as functions of the density. The 
hydrodynamic (1), intermediate (2), and ideal-gas (3) ranges 
are shown, (e) One-body density matrix pi(r) at n = 1 and 
large r: DMC calculation (dashed curve), hydrodynamic cal- 
culation (|17p (solid curve), and asymptotic form (|18|) for the 
infinite system (dot-and-dash curve). The endpoint r/^/n = 5 
of the abscissa axis corresponds to the boundary L/2 of the 
system. 



to the MC data. The results are described well (within 
the limits of 0.005) by the following polynomial fitting 
curve: 

n Q /n = a l exp(& l lnn + Cq In 2 n + dg\n 3 n), (19) 



where = 0.3822 

d n - 
a - 



h n - 

°o — 



0.2342, Co 1 = -0.02852, and 
—0.001594. According to the data presented in 
Fig. 3(c), no/n « 1/3 at n = 2, and, at n = 8, the frac- 
tion of the Bose condensate amounts to only no/n w 1/5. 
This is indicative of strong correlations in the dense gas of 
2D dipolar excitons in CQWs in their experimentally ob- 
served collective state. It should be specially noted that, 
at so small fractions of the condensate, the application of 
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the mean-field methods (the Gross-Pitaevskii mean-field 
approach), the perturbation theory calculations, as well 
as the Bogoliubov approximation, is unjustified and can 
be used only for qualitative purposes. 

We also note that, in the macroscopic limit (N — > oo), 
the Bose condensate fraction is even smaller than that at 
N = 100. Thus, at n = 2 and n = 8, respectively, we 
found that uq/ti » 0.29 and no/n ~ 0.17. 

Figure 3(b) shows that hydrodynamic expression (fT7|) 
excellently coincides with the MC calculation of p\{r) 
at distances r longer than the characteristic microscopic 
phonon length scale ro (i.e., at r > ro). However, at 
r < 7*o, this coincidence vanishes very rapidly. Therefore, 
the crossover range of the domain of applicability of hy- 
drodynamics for 2D dipolar excitons turns out to be very 
narrow. 

The behavior of the characteristic linear microscopic 
phonon scale ro of the system in relation to the density 
is shown in Fig. 3(d). It is clearly seen that, in the en- 
tire density range, this scale has the order of the average 
interexciton distance, 



1 



ro 



(20) 



In this case, the healing length £ = 1/^/2/ii proves 
to be considerably smaller than ro at all densities (see 
Fig. 3(d)). This contradicts the weak correlational be- 
havior of 2D dipolar excitons. 
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Figure 4: (a) Exciton pair distribution function at T — for 
different densities n = 2'~ 9 (i = 1,2. ..12). The inset shows 
density dependences for the diameter a of the dipole effective 
hard disk, polynomial fitting curve (|22f) . energy dependent 
dipoie scattering iength, a s (|23|) . and cfassical turning point 
for the dipoie potential, a cl = l/^/2E /N. (b) VMC and 
DMC results and their iinear extrapolation (LE) for pair dis- 
tribution function at n = 1 (curves practically coincide) . The 
arrow clarifies the definition of the diameter a of the dipoie 
effective hard disk. 



C. The pair distribution. The diameter of the 
dipoie effective hard disk. The influence of the 
exciton internal structure. The energy-dependent 
scattering length 

The polar-circle-averaged pair distribution function of 
excitons, 



P2 (r) = f (p(r)p(0))d^/2n, 
Jo 



(21) 



where p(r) = ^> + (r)^(r) is the exciton density operator, 
simulated at different densities n is shown in Fig. 4(a). 

At high densities (n > 1/8) , the pair distribution func- 
tion exhibits a clearly pronounced hump, corresponding 
to a short-range order. At very high densities, there are 
also weaker humps. This is indicative of strong corre- 
lations in 2D dipolar exciton system in CQWs in the 
low-temperature exciton phase studied in [H, [y, 0, @, 0] • 

The inset of Fig. 4(a) shows the density dependence of 
the diameter a of the effective hard disk of a dipoie. We 
define this diameter as the distance from the origin to the 
point of intersection of the short-wavelength tangent line 
to the pair distribution function with the abscissa axis 
(see Fig. 4(b)). The calculated points are described well 
(within 2%) by the following polynomial fitting curve: 

a = a a exp(6 a Inn + c a In 2 n + d a In 3 n), (22) 



where a a = 0.293, b a = -0.324, c a = -0.0123 and d a = 
0.00096. 

The calculation of the diameter of the effective hard 
disk makes it possible to evaluate the influence of the in- 
ternal exciton structure, which is not taken into account 
in the model of dipolar excitons. If the exciton diame- 
ter is smaller than the diameter a of the dipoie effective 
hard disk, the neglect of the internal structure of dipolar 
excitons is justified. This is connected with the fact that 
the probability density of one exciton to be at a distance 
of r < a from another exciton, which is proportional to 
P2O"), is small (Fig. 4(b)). 

The average electron-hole distance a e h in the exciton 
can be taken to estimate the exciton diameter in GaAs. 
In structure (1), which consists of GaAs CQWs [1, 0], 
the average distance a e h = 17 nm; therefore, for the di- 
mensional exciton density n < 3.2 • 10 10 cm -2 , a e h < 0,, 

If the exciton diameter is greater than a, the internal 
structure somewhat affects the microscopic properties of 
dipolar excitons. However, up to densities n correspond- 
ing to the boundary between the strong and weak cou- 
pling regimes, where the dipoie barrier ceases to suppress 
exchange effects (see Sections HI and HIj) . this influence still 
can be neglected. 

In structure (1) [B|, 0], at the maximal exciton den- 
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sity 5.17 • 10 10 cm -2 (see Eq. ([6])), corresponding to the 
boundary between the strong and weak coupling regimes, 
we find a = 14.2 nm w 0.840^. 

Finally, in the inset of Fig. 4(a), we present the scat- 
tering length a s = a s (p), 



a s (p) = a a s exp{K hip + c a s In 2 p + d a s In 3 p) 



(23) 



of two dipoles with the relative momentum p ~ y/ 2E / N 
and the energy of order of Eq/N. Here, a a s = 0.68845, 
b a s = -0.45897, c a s = -0.03098, and d a s = 0.002096. Fit 
([23]) is based on the calculation of Ref. [59| and, at 
0.1 < p < 20 (i.e., at 1/340 < n < 9.3; see Eq. ((Til ), 
coincides with it within 0.04%. It is seen from the inset 
that, at low densities, the diameter a of the effective hard 
disk is smaller than the scattering length a s . 



D. The static structure factor. The sound velocity 

The structure factor of 2D dipolar excitons is defined 
via the Fourier transform of their pair distribution func- 
tion p 2 (r) = (p(r)p(O)) as 



5(p) = 1 



ipr 



(p 2 (r)/n- n)dv. 



(24) 



For a system with a finite size L and with periodic bound- 
ary conditions, this parameter is a discrete function of 
the momentum (p = (2ir/L)n, where n = {n x , n y } is the 
integer-valued vector). However, the experiment is, as a 
rule, interested in large 2D Bose condensed exciton sys- 
tems with a rather large exciton number and a complex 
boundary. For their description, it is convenient to use 
the expression for the structure factor in the correspond- 
ing infinite system. 

For momenta p 3> 2n/L, the form for the structure 
factor of the infinite system, S(p) = S(p), can approxi- 
mately be obtained from Eq. ([24]) for a flat trap of the size 
L (with a number of excitons equal to, e.g., N — 100), if 
Eq. (I24p is smoothed by integrating over the polar angle, 

S(p) = 1 + J J Q (pr)(p 2 (r)/n - n)2iTrdr, p > 2n/L. 

(25) 

Here, it is taken into account that J Q e~ lxcosv dip — 
2irJo(x), where Jq(x) is the zeroth-order Bessel func- 
tion. For a density of, e.g., n = 1, the difference between 
smoothed structure factors (l25|) determined for two val- 
ues of the exciton number, N = 100 and N = 200, does 
not exceed 0.35%. Consequently, smoothed form (|25|l for 
the structure factor approximates well this parameter for 
the infinite system at p 3> 2ir/L. 

However, at very small momenta, p < 2nh/L, at which 
effects of the finiteness of the system manifest themselves, 
smoothing (|2"5")1 does not yield the correct result for infi- 
nite systems. To obtain the correct result for an infinite 
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Figure 5: (a) Exciton structure factor at T — for different 
densities n = 2 1-9 (i = 1,2. ..12). (b) Sound velocity deter- 
mined from the slope of the structure factor (squares) and 
based on fitting (|11[1 to the energy (solid curve) as a function 
of the density. 



system at such small momenta, it is necessary to join 
Eq. (|25[) with the long- wavelength hydrodynamic asymp- 
totics for the structure factor known from the general 
theory. For higher accuracy, we will use the third-order 
long- wavelength asymptotics in momentum f&oj ] . 



S(p) = sp(l + 7p 2 ), p < 2tt/£, 



(26) 



valid on hydrodynamic scales p <C Po = 1 /tq ~ ^fn (see 
Section IIII Bl) . Here, s — const > is the parame- 
ter responsible for the long-wavelength behavior of the 
structure factor, while the quantity 7 = const deter- 
mines the third order in p. The smooth joining (the 
coincidence of the functions and of their first deriva- 
tive) of the smoothed structure factor (|2"5"|) and the long- 
wavelength asymptotics (|26[) at intermediate momenta 
2n / L <C p <C Po turns out to make it possible to rather 
precisely determine the value of s. For N = 100, the 
prediction for s, s = l/^/Snn^, based on the polynomial 
fitting (HI]) to the energy (see Eqs. l[9>([Ti|). (|27|l. (|2~gjl) 
yields the difference from the macroscopic limit that is 
less than by 0.08%. This is indicative of the adequacy of 
our approximations. 

Figure 5(a) depicts the structure factor of the cor- 
responding infinite system calculated by Eqs. ([23)1 and 
(f2"6")> for different densities. At high exciton densities 
(n > 1/8), the curves exhibit a hump corresponding to a 
short-range order. At densities 1/4 < n, < 4, this hump 
is clearly seen. This testifies to strong correlations in 
the dense gas of 2D dipolar excitons in CQWs studied in 

Figure 5(b) shows the dependence of the sound velocity 
c s in the collisionlcss regime on the density, which we 
determine from hydrodynamic considerations from the 
parameter s, 



1 

2s 



(27) 



(see Eq. (|29l) below). The calculated points for c s fall 
well on the curve corresponding to the expression of this 
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quantity via the adiabatic compressibility, 



i/x = \/2Tm( 



(28) 



(the dimensionless adiabatic compressibility C is deter- 
mined by Eqs. ©, ([TUJ) , and dH|) from fitting fll]) to the 
energy; see Section fill A[) . This indicates that our DMC 
simulation is internally consistent. 



E. The excitation spectrum. Characteristic 
microscopic scales of the system 

In the collisionless regime, the gas phase of the Bosc 
system is known to have only one spectral branch. There- 
fore, the spectrum of elementary excitations of Bose con- 
densed 2D dipolar excitons in an infinite system at T = 
can approximately be calculated by the Feynman formula 



p 



2S(p) 



(29) 



which is valid in the long-wavelength (hydrodynamic) 
range (r ro, p <C po), where the excitation spectrum 
corresponds to phonons. Approximately, we have (see 
Section lmDj) S(p) « sp(l + jp 2 ) and e p « c s p(l - -fp 2 ). 
Formula (j2"9")l is also true in the short-wavelength range 
corresponding to an ideal gas (r <C nc p ^ Pig = 
1/ric), where S(p) « 1 and e p w p 2 /2. However, on 
intermediate scales (rjc <C r <C ro, pjg >p> Po)> the 
Feynman formula (f2"9"| is quantitatively not valid. Never- 
theless, it still can be used for qualitative estimates. The 
details are discussed in Appendix B. 

In Figs. 6(a) and 6(b), the excitation spectra in an in- 
finite gas of 2D dipolar excitons analytically calculated 
by Eq. (f2"9"|) for different densities are shown. It is seen 
that, in the entire density range used in the calculations, 
the excitation spectrum is far from the Bogolyubov shape 
Ep = -\/p 4 /4 + c 2 p 2 (although, at low densities, it qual- 
itatively resembles the Bogolyubov spectrum). At high 
densities (1/4 < n < 4), effects of strong correlations are 
clearly seen. At the maximal density used in the calcula- 
tion, n = 8, a roton minimum appears in the spectrum. 

The ranges of long- wavelength (hydrodynamic) , short- 
wavelength (ideal-gas), intermediate, and microscopic 
scales are clearly seen in Fig. 6(c). The diagrams of the 
long- wavelength, intermediate, and microscopic ranges 
were presented above in Fig. 3(d) (see Section UlI Bp . 



F. Temperature dependence of the local superfluid 
density. Quasicondensation and BKT transition 
temperatures 

In the quasicondensed [62| phase of excitons far from 
the crossover between the quasi-classical regime and the 
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Figure 6: (a) Excitation spectrum at T = for different densi- 
ties n = 2 1-9 (i = 1, 2. ..12). (b) Comparison of the excitation 
spectra with the Bogolyubov shape, (c) Excitation spectrum 
for the density n = 4 (solid curve) and the spectrum of the 
ideal gas (dashed curve). The hydrodynamic (HD; r ^> ro, 
p <C Pa), ideal-gas (IG; r ria, P 2> Pig), intermediate 
(ro > r » riG, Po < P < Pig), and microscopic (r < r , 
p 2> pa) ranges are shown. 



quasicondensed phase, elementary excitations form a 
nearly ideal gas. In this case, the temperature depen- 
dence of the fraction of the local (vortex-unrenormalizcd 
[63j ]) superfluid component of quasicondensed 2D dipolar 
excitons can approximately be calculated by the Landau 
formula 16411 . 



ni(T) 



p 2 dn(Sp) pdp 
T del 27rn9eXl 



T q -T~T q . (30) 



Here, n; is the local superfluid exciton density (in the g e 
spin degrees), 



n(ej) 



1 



- l 



(31) 



is the Bose distribution of the ideal gas of elementary 
excitations with the spectrum at the exciton temper- 
ature T, and T q is the estimate for the temperature of 
quasicondensation crossover, at which the local super- 
fluid component and the local long-range order [57] van- 
ish gradually. 

The global superfluid component (which takes into ac- 
count the vortex renormalization [24, 63]) was calculated 
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ni{T)/n = 1 - (1 - ni{T c )/n){T/T c y>, (33) 



7; = + bl In n + c7 In 2 



+ dj In 3 n, (34) 



where a] = 4.186, 6/ = 0.5825, c] = 0.0941, and 
dj = 0.00706, while T c < T' is the temperature of 
the BKT superfluid transition [23, 24 , 25|, at which the 
global superfluid component in the infinite system van- 
ishes jump- wise 0, For T c , we use the equation 
( (see also Appendix B) 



24 
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in which the dielectric permittivity e of vortex pairs at 
the BKT transition can be taken from the 2D x-y model 
(e = 1.1349 ± 0.0005; see 0). 

In Figs. 7(b)-(d), we present the analytically calcu- 
lated points (denoted by squares) for the temperatures 
of quasicondensation, T q (see Appendix B), and of the 
superfluid BKT transition, T c , as well as for the local su- 
perfluid fraction, ni(T c )/n, at the BKT transition in the 
infinite system. These points fall very well on the cor- 
responding polynomial fitting curves (within the limits 
0.0025, 0.0005, and 0.0025, respectively), 



Figure 7: (a) Temperature dependence of the local super- 
fluid exciton fraction in the g ex — 4 spin degrees for dif- 
ferent densities n = 2'~ 9 (i = 1,2. ..12) (solid curves), low- 
temperature power-law fitting curves determined by Eq. (|33[) 
(dashed curves); and extrapolation to the range of the qua- 
sicondensed crossover (dotted curves), (b-d) Density depen- 
dences (squares) for (b) the quasicondensation temperature 



T q /T 



-b q lnn- 



-Cq" In 2 n- 



■d q ln J 



-e q In n, 



T c /T Q = a T c +b T c In) 



T i 2 

s„ m 



di In 



n, 



ni (T c ) /n = a c + b c In n + c[ 



In 2 n 



In 3 n. 



0.038, c 

-0.001146, and e 1 ' = -0.000086: 



T q , (c) the BKT transition temperature T c , and (d) the local Here a T = 5344 
superfluid fraction at T — T c ; solid and dashed curves corre- 
spond to the polynomial fitting and to phonon-model calcula- 
tion, respectively, (e) Exponent 7; in Eq. (|33|l as a function of 
the density: (squares) the fit for each density, (solid curve) the 
general fitting by Eq. (|34[) . and (dashed curve) the calculation 
within the phonon model. 



-0.00395. 
= 0.2146, 



(36) 
(37) 
(38) 



4 = 

b T c = 
1 _ 



m 



44j | . To calculate the local component, one can set 
„ ~ e p in Eqs. (001) and §^ . The details of 

the calculation of the local superfluid component are dis- 
cussed in Appendix B. 

The dependence ni{T)/n analytically calculated by 
means of Eqs. (f2"9"|) - (j3"Tj) with = e p at different densi- 
ties is shown in Fig. 7(a), where the quantity 

T = g ex Tdeg = 27m (g ex = 4 in GaAs). 



(32) 



0.00414, c T c = -0.000877, and <F C = -0.0000446; a u 
0.9745, b\ = 0.0188, c\ = -0.00398, and d' c = -0.000202. 

According to our results, the phonon model, which is 
widespread in the literature [5^, [3(| (in which e p = c s p) , 
describes rather well the temperature of the BKT transi- 
tion (Fig. 7(c)) but fails to adequately describe the tem- 
perature of the quasicondensation crossover (Fig. 7(b)). 
(The phonon model does not take into account the ro- 
ton bend of the spectra on intermediate scales (see 
Figs. 6(a),(c)), which makes the main contribution at 
high densities and T ~ T q .) 

It is important that, at high densities (1/4 < n < 
4), the local superfluid fraction at the BKT transition, 
ni(T c )/n, is close to unity. This is inductive of strong 
correlations in the dense 2D dipolar exciton gas in CQWs 



denotes the degeneracy temperature for the correspond- 
ing spin-polarized excitons. 

In the temperature range < T < 1.2T C , with an error 
less than 0.003, this calculation of the local superfluid 
fraction coincides with the following power-law fitting 



correlations m tne dense zu amc 
[1 00,0,11 (see also Ref. Q). 



In addition, for the spin-depolarized dense 2D dipo- 
lar exciton gas in GaAs (g ex — 4), the temperature of 
quasicondensation T q ~ 0.5T (see Fig. 7(b)) is approxi- 
mately two times higher than the temperature of degen- 
eracy Tdeg = T /g e x = 0.25T (see Eq. (|3"2l ). whereas the 



11 



BKT transition temperature T c w 0.2T (see Fig. 7(c)) is 
only slightly lower than T ( i eg - This strongly differs from 
models widespread in the literature of ideal and weakly 
correlated 2D exciton gases in QWs for which these tem- 
peratures are logarithmically small compared to the de- 
generacy temperature T^ eg [53| . 

Therefore, the condensed state of 2D dipolar exci- 
tons in CQWs and SQWs can be experimentally obtained 
much simpler than in the case of weakly correlated exci- 
tons. 



IV. A HARMONIC TRAP AT T = 

Modern experimental methods are capable of harmonic 
trapping of 2D dipolar excitons in the plane of a quantum 
well. This trapping can be obtained with the help of an 
inhomogeneous compression of the sample, caused by the 
pressure of a tip on its surface ll|, LL2| , as well as with 
the help of an inhomogeneous electric field in electrostatic 

traps d, e, Ei, m ei m . 

In the latter experiments on Bose condensation of ex- 
citons, their number in a harmonic trap is large, N ^> 1 
|1, [HI, [13, EH- In this case, the phonon microscopic 
1/V^(0) (see Eq. O; n(R) is the 



length scale ro 
total exciton profile in the trap) proves to be signif- 
icantly smaller than the trap size Ltf ~ \/N/n(0) 
[L TF = yj2N/im{Q) is the Thomas- Fermi (TF) radius). 
In the calculated density range, other microscopic scales 
(£, tig, a, etc.) are even smaller than ro (see Section 
IIIII) . Therefore, microscopic properties of excitons in a 
trap can be formed on local scales, along which the po- 
tential of the trap varies continuously. 

Moreover, if the exciton number N in a trap is substan- 
tially greater than Air 2 w 40 the hydrodynamic range 
p <C 1/ro ~ y n(0) considerably exceeds the momen- 
tum discreteness step 2tt/Ltf ~ 2iry/n(0)/N <C yn(0), 
which is determined by boundary effects in the homoge- 
neous system with the size L = Ltf (see Sections MI Bl 
and IIII D|) . Therefore, in the exciton system in a large 
trap containing N ^S> 40 excitons, the sound range of mo- 
menta, 2-7T / 'Ltf <C p <C 1/fo does exist and is not masked 
by finite-size effects, arising on scales p < 2tt/Ltf- 

Therefore, a sufficiently large harmonically trapped 2D 
dipolar exciton system at T = can be fairly well de- 
scribed in the local density approximation (LDA) [69J] . In 
this approximation, the local microscopic exciton prop- 
erties near the point R of the trap are approximately 
replaced by the properties of the corresponding homo- 
geneous system with the density n equal to the exciton 
density n(R) in the trap at the point R. 

Thus, within the framework of the LDA, the total ex- 
citon density profile in a (symmetric) trap n(R) is deter- 
mined from the TF equation, which, at T = 0, takes the 
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Figure 8: (a) Total (n(R)/n(0)) and Bose condensate 
(no(r)/n(0)) exciton profiles in a harmonic trap in the LDA at 
T = and Q' e = 2e 2 D/e = 10 for different densities at the cen- 
ter (n(0) = 2 4 " 9 , i = 3, 4.. .12). (b) Comparison of the Bose 
condensate profiles at different n(0) with the Thomas-Fermi 
inverted parabola. 



form 



V(R) + n(n{R)) = n, 



(39) 



where V(R) = (m/2)uj 2 lo R 2 is the potential of the trap 
with the oscillator frequency Who, M is the chemical po- 
tential of excitons in the trap counted from the exciton 
band edge, and fj,(n(R)) is the local chemical potential of 
excitons in the trap equal to their chemical potential in 
the homogeneous system with the density n = n(R). 

The density of exciton Bose condensate at T = in a 
harmonic trap can be calculated in the LDA as 



n (R) = n Q (n(R)), 



(40) 



where no(n(R)) is the Bose condensate density in the 
homogeneous system with the total density n = n(R). 

Figure 8(a) shows the total and Bose condensate ex- 
citon profiles in a harmonic trap at T — analytically 
calculated in the LDA for different densities n(0) at the 
trap center. The quantities n{n) and no(n) in the infi- 
nite homogeneous system were approximately calculated 
using the corresponding fits (jTTJ) and (fT9"|) (see also ([5]) 
and ([5])) for the homogeneous finite system consisting of 
N = 100 excitons. 

It is clearly seen from this figure that, at C' e = Me/^o = 
2e 2 D/e = 10, the total density n(R) for all n(0) nearly 
ideally fall on the profile of the inverted TF parabola (see 
Fig. 8(a)), 

n(R) « n(0)(l - R 2 / L 2 tf )6{L T f - R) (41) 
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where 9(x) = at x < and 9{x) = 1 at x > 0. How- 
ever, this fact is determined not so much by a weak non- 
linearity of the collisional contribution to the chemical 
potential as by a large value of the electrostatic contri- 
bution /i e = CgT = lOXo, which is linear in the density 
(see Eq. (I3"2"10 . (The values of the electrostatic contribu- 
tion to the chemical potential CL ~ 10 were realized in 
experiments @, % 0, i, @, H [13,1131 •) 

However, the shape of the Bose condensate profile 
no(R) at high densities (1/4 < n < 4) appreciably 
differs from that of the inverted parabola no(0)(l — 
R 2 IL 2 TF )6{L T F - R), Fig. 8(b). This is evidence in fa- 
vor of strong correlations in the dense 2D dipolar exciton 
gas. 

Note that the LDA makes it possible to predict the 
frequencies of collective oscillations of the compression 
mode, when the excitation is caused by a sharp change 
in the frequency of the trapping potential [7(j. The fre- 
quency of this mode depends on the particular form of the 
two-particle interaction (in contrast to the mode that is 
caused by the displacement of the center of mass and that 
depends only on the trap frequency). This makes it pos- 
sible to experimentally investigate the equation of state 
(see Eq. (I11|) V Another important quantity measured in 
experiment is the release energy. A rapid switching-off 
of a trap turns the trapping potential to zero but practi- 
cally has no effect on the kinetic energy and the part of 
the potential energy related to the pair interaction. The 
released energy can be measured during the spread. Us- 
ing the LDA. it is possible to estimate the release energy 



V. THE POSSIBILITY OF THE 
EXPERIMENTAL OBSERVATION OF STRONG 
EXCITON CORRELATIONS 

Effects of strong correlations in a dense 2D dipolar 
exciton gass in CQWs can be found upon observation 
of the particular features of exciton luminescence in a 
longitudinal magnetic field. 

Indeed, in the absence of a magnetic field, according to 
the momentum conservation for one-photon exciton re- 
combination, the momentum p of a recombining exciton 
is equal to the projection of the momentum of an emit- 
ted photon onto the QW plane, (Huj/co) sm9. Here, 9 is 
the angle between the emitted photon and a normal in a 
vacuum, cq is the light velocity in a vacuum, and u> is the 
photon frequency (now we pass to the ordinary, dimen- 
sional, units). However, if a longitudinal magnetic field 
H\\ is applied to the system, the dispersion curve e p of 2D 
dipolar excitons is shiftedby the quantity pu = &DH\\/cq 
[71| (this effect is equivalent to the excitation of diamag- 
netic currents in a system of coupling e and h by a longi- 
tudinal magnetic field that was described in [l[). As 
a result, the relation between the emitted-photon angle 9 



and the recombined-exciton momentum p takes the form 



(Huj/cq) sin( 



PH • 



(42) 



If we are interested in luminescence along the normal 
(9 = 0), then, p = pn (see Eq. (|l2"]l). so, e p = e PH . 
Hence, for the spectral-angular luminescence along the 
normal in the field H\\ at T = we obtain 

lHo(Lu)=lii 5(Lu-n + e PH /h), ph^>Pt, (43) 

where IgL is the spectrally integrated angular lumines- 
cence along the normal in the field Hu , Q is the exciton 
resonance frequency (for the luminescence at Hu = 
and 9 — 0), and px — T/c s is the characteristic ther- 
mal momentum of Bose condensed excitons, which sepa- 
rates the ranges of thermal p <C pt and zero-temperature 
Pt "C p -C Po quasicondensate phase fluctuations • 

Eq. (l43|) ) shows that the measurement of the lumines- 
cence line redshift [73| e PH /h in the magnetic field H\\ = 
coPh/zD makes it possible to directly determine the ex- 
citations spectrum e p . At D = 12.3 nm, m — 0.22mo, 



12.5 



2-10 10 cm" 2 , and H 



T, we have 



p H /h = 1-49 • 10 6 cm" 1 = 10.5Vn w p IG /h 

(see Section UT1 and Fig. 3(d)). Therefore, in structure 
(2) 0, Q containing an exciton gas with the density n — 
2 • 10 10 cm -2 , the field H\\ = 8 T covers both the sound 
(hydrodynamic, p -C p Q ) and the intermediate (po "C 
P *C Pig) spectral ranges (see Fig. 6(c)). In this case, 
at T = 0.2T = 1 K, m = 0.22m , and C = 3.7 (see 
Fig. 2(b) and Sections HT1 IIII A| ) . the thermal momentum 
Pt corresponds to the field 



H 



eD 



(To 



0.2 T. 



Knowing the slope of the measured spectrum e PH , one 
can calculate the sound velocity c s f=a £ PH /ph (ph *C po; 
see Fig. 6(a)). This yields the dimensionless adiabatic 
compressibility (see (j28|l ) 

C = m 2 c 2 s /2TTh 2 n. 

Alternatively, the adiabatic compressibility £ can be 
found from the measurement of the spectrally integrated 
luminescence along the normal in the magnetic field, 

1 e=o- 

Indeed, the spectrally integrated luminescence into 
the solid angle dtt = sin 9dipd9 near the normal is 
connected with the momentum distribution of excitons 
n PH = J e~~fi pHr pi(r)dr (i.e., the number of excitons 
with the momentum pu ) by the following relation: 



^d^l — hi 



— ) 

TO 



"PH 



dpn 

(2nh) 2 



(44) 



Here, To is the lifetime of an isolated spin-depolarized 
exciton in the ground state and k = 1/2 (k = 1) in the 
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case of the measurement of the luminescence only on one 
side (on both sides) of the QW plane. 

On scales pt <C ph <C Po that are the most interesting 
to us, the momentum distribution at T = is given by 
the Fourier transform of Eq. CR 



no mc s 
n 2p H ' 



Pt *C ph <C po- 



(45) 



From Eqs. (|44|) and (|45|) . we find the sought relation be- 
tween the angular luminescence along the normal, I^L , 
and the dimcnsionlcss adiabatic compressibility £, 



2 e=o 



Io 



qf c o 



4ireD\/mTo H\\ ' 



Pt < eDH\\/co < Tia/™- 
(46) 

Here, Jo = K(Ml/ro)no is the exciton Bose condensate 
luminescence at T = 0. At n = 2 - 10 10 cm" 2 , L> = 12.3 
nm, T = 0.1T = 0.5 K, m = 0.22m , e = 12.5, and 
C = 3.7, the momentum range pt <C eDHu/co <C Po m 
Eq. (j4"B"|) corresponds to the magnetic-field range 



cpT 
eD 



(T 



< Hu < 



or 0.1 T < < 0.75 T. 

The luminescence intensity of the Bose condensate at 
T = appearing in Eq. (I46|) can be approximately cal- 
culated as the contribution of thermal phase fluctuations 
to the quasicondensate luminescence at low temperatures 
(see Eq. dSl), 



H 2np H dp H 



Qr 



m 

TO 



T « T g , (47) 



where q r = Ml/co is the radiation bandwidth in a vac- 
uum. Here, it is taken into account that the exciton gap 
Kl ~ 1 eV considerably exceeds all characteristic energy 
scales for excitons s PH ,T < 1 meV <C fiil, so, Tilo « M7. 
It is also taken into account that dpn — q r dsiri9 = q r d6 
at p H = p and 9 = (see Eq. (JUJ). (Eq. (|47]) can be ob- 
tained by extending the method described in Appendix 
A to the case T ^ 0, see Refs. [H 

Finally, by measuring the quantity (|4"?| . one can find 
the important parameter of exciton correlations, namely, 
the Bose condensate density at T = 0, 



n 



ft 



VI. CONCLUSION 

By means of the ab initio simulation and analytical 
calculations, we studied in detail the microscopic proper- 
ties of Bose condensed superfluid 2D dipolar excitons in 
QWs. For a homogeneous exciton system, we numerically 
calculated the ground-state energy, the one-body den- 
sity matrix, and the pair distribution function. Based on 



these numerical calculations for the homogeneous system, 
we analytically found the structure factor, the excitation 
spectrum, the Bose condensate density, the temperature 
dependence of the local superfluid density, the chemical 
potential, the adiabatic compressibility, the sound veloc- 
ity, and the characteristic microscopic length scales (the 
phonon scale, the healing length, the scattering length, 
the diameter of the effective hard disk of a dipole, etc.), as 
well as the quasicondensation crossover and Berezinskii- 
Kosterlitz-Thouless transition temperatures. For a har- 
monic trap with a large excton number, we analytically 
calculated the total and Bose condensate profiles. 

We showed that, in all experiments on Bose condensa- 
tion performed at present, 2D dipolar excitons in coupled 
QWs prove to be strongly correlated. As a consequence, 
exciton Bose condensation and superfluidity can experi- 
mentally be achieved much easily than this is predicted 
by the models of an ideal and weakly correlated exciton 
gases. 

The results obtained for the ground-state energy 
Eq/N, chemical potential a = fa + fx e , the adiabatic 
compressibility %, the sound velocity c s , the excitation 
spectrum e p , and the Bose condensate no and local su- 
perfluid n/(T) densities, as well as for the characteristic 
microscopic phonon length scale ro , can be used as refer- 
ence data in construction of quantitative hydrodynamic 
theory of the optical properties of 2D dipolar excitons 
in QWs, in particular, of their luminescence, coherence, 
and nonlinear effects. 



Appendix A: LONG- WAVELENGTH 
HYDRODYNAMIC ASYMPTOTICS FOR THE 
ONE-BODY DENSITY MATRIX 

On long- wavelength scales (r ^> ro), the one-body den- 
sity matrix (| 16[) of Bose condensed 2D dipolar excitons 
at T = can be calculated within the framework of the 
hydrodynamic method in quantum field theory [22| , 



Pi 



(r) = Ce _<(l?w_#(0))2/2> = 



n e 



(0(r)0(O)> 



r > r . 
(Al) 

Here, <p(r) is the phase of the exciton field operator ^(r) 



(\l/(r) = e 1 ^^ y/ /5(r)), which satisfies the commutation 
relation [<p(r),(p(r')] = 0, and C — n e^ < -°^ 2 ' > = const 
does not depend on r. We assume that an ultraviolet cut- 
off is imposed (in this case, the quantities C and ((<p(0)) 2 ) 
diverge in the ultraviolet limit, but, at r 3> ro, no and 
(tp(r)tp(O)) are finite in the large-momentum limit). 

At T = 0, the phase-phase correlator (<£(r)<£(0)) in 
Eq. (|Aip can be expressed via the phase-phase Green's 
function (in dimensional units), 



D^(w) = -i /<T[0(r,t)0(O, ())]>< 
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as follows: 

(^(0)> = hE 



iD^{oj)^. (A2) 

Z7T 



where T is the chronological operator. 

In the hydrodynamic (long-wavelength) approxima- 
tion, in which the hydrodynamic Hamiltonian acquires 
the form local in time, the function D^ v (u)) is given by 



Hme 2 /n s p 2 
'?uj 2 — E 2 + id ' 



(A3) 



where n s is the global superfluid exciton density (in the 
g ex spin degrees). 

In the calculated density range {nx^ < 8), excitons 
are in the gas phase [36j]. In addition, we neglect the 
external random potential, which is unavoidably present 
in heterostructures. In this case, the global superfluid 
exciton density at T = coincides with the total density 
(see [H| and Section UTTF)) . 



n s = n. (A4) 

From Eqs. (|A2|) - J A4|) we obtain the following equation 
for the correlator (<^(r)<^(0)) at T = 



<^(0)) = ^E p#o eiPr 



me p 
2np 2 



-pr 



me 



p 



1 



V e iP' 

p 



e *Pr ( HlfP- 

^ \ 2Np 2 

p#0 v F 



r > r . 



1 

AN 



1 



(A5) 



Here, we used the fact that, at r > r (i.e., at r ^ 0), 
the 5-like contribution 



AN ^ 



-E 

An ^ 



S(r - nL) 



is zero in a box of a size L (n = {n x , n y } is the integer- 
valued vector). 

By substituting (|A5|) into (|A1[) and averaging over the 
polar angle of the vector r, we arrive at the sought hydro- 
dynamic equation (|17[) for the one-body density matrix 
of Bose condensed 2D dipolar excitons at T = 0, which, 
on long- wavelength scales (r ^> r ), coincides excellently 
with our DMC simulation (see Fig. 3(e)). 

In conclusion, we note that, in interacting systems 
with the finite scattering length, the spectrum of excita- 
tions £ p , at p — > co, exponentially approaches the spec- 
trum of the ideal gas e p = p 2 /2m + 0(e~ p / pi ) (where 



pi = const > 0; see also item (2) in Appendix B). Con- 
sequently, at p — > oo, the quantity me p /2Np 2 — 1/AN 
under the summation sign in (|A5[) is of the order of 
0(e _p / pi ). Therefore, Eq. (|A5p for the phase-phase cor- 
relator converges in the ultraviolet limit. This means that 
the well-known ultraviolet cutoff [62l | in quantum-field 
hydrodynamics at the boundary of the long-wavelength 
band (r ~ ro, p ~ pa) proves to be redundant in (|A5[) : 
an infinitely small cutoff suffices to ensure convergence. 

However, one should bear in mind that the hydro- 
dynamic (long-wavelength) approximation fails to ade- 
quately describe the phase-phase correlator on micro- 
scopic scales r -C r Q , p ^> p . So, the range p 3> Po 
is incorrectly taken into account in (|A5|) . As a result, the 
hydrodynamic approach used in this study introduces a 
certain error in the calculation of the one-body density 
matrix and, therefore, in the Bose condensate density. 
However, direct calculation shows that, if a continuous 
finite ultraviolet cutoff is introduced in (|A5[) at p ~ po, 
the correction to the Bose condensate density will rapidly 
vanish at large N. The numerical calculation of the sum 
in (1A5|) shows that, at N = 100, this correction is of the 
order of 0.001. 



Appendix B: APPLICABILITY CONDITIONS 
AND CALCULATION DETAILS 

Here, we will discuss subtle aspects of calculations of 
the excitation spectrum, the local superfluid density, and 
the quasicondensed and superfluid transition tempera- 
tures. 

(1) As Boronat et al. showed [6l|, in strongly corre- 
lated systems on intermediate scales, the true spectrum 
lies lower than the spectrum calculated by the Feynman 
formula (although, the weaker the correlations, the more 
exact the Feynman formula). Therefore, in the strong- 
correlation regime, our calculation of the local superfluid 
fraction, as well as the quasicondensation T q and the 
BKT transition T c temperatures is valid only qualita- 
tively. 

(2) For simplicity, the excitation spectrum was esti- 
mated omitting the effects connected with a large damp- 
ing of elementary excitations above the endpoint of the 
spectrum 

(3) The Landau formula is valid only in the collision- 
less (not hydrodynamic) regime when elementary excita- 
tions weakly interact, and their damping is small. This 
is definitely not the case above the spectrum endpoint 
p c . In a dense system, the spectrum endpoint lies high 
(see Fig. 6(a) and (Si), so that the range p > p c is cut in 
(|30|) by the Bose exponential ([31]) . At low densities, the 
spectrum endpoint p c corresponds to smaller frequencies, 
but a significant damping appears only above the sound 
range of the spectrum. This range is also cut by the Bose 
exponential (see (|3"0"|) - (l3"2"|) and Figs. 6(a), 7(b)). However, 
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at very low densities n, where excitons are weakly cor- 
related, the hydrodynamic regime can occur precisely in 
that range of the spectrum which makes the main contri- 
bution to the Landau formula. In this case, Eq. (|30[) docs 
not hold. However, according to our direct estimation, 
even at a smallest calculated density of n = 1/256, such 
a situation is not realized. 

(4) In calculation of the local component, we set ej » 
£p =0 — £ p in Eqs. ([50)1 and ([51]) . Indeed, at low tem- 
peratures, T <C T q , the Bose exponential ([51)1 leaves in 
(|30| only the contribution from long- wavelength scales, 
on which w c s (T)p w c s (0)p w e p . At high tempera- 
tures, (T - T„ T < T g ), owing to the factor (p 2 /2)pdp 
in Eq. lp3Qj) . long- wavelength scales do not contribute. 
On scales in the range of shorter (microscopic) wave- 
lengths, the spectrum calculated by the Feynman formula 
is weakly affected by temperature variations. The latter 
fact was revealed in the course of the Monte Carlo simu- 
lation of the pair distribution function and the structure 



factor in liquid helium up to temperatures so high (4.2 
K) as the doubled temperature of the A point, T\ = 2.1 
K (z|. 

(5) The temperature of quasicondensation T q we for- 
mally calculate from the condition ni(T q )/n = 0, where 
ni(T)/n is given by Eq. ([50)1 with = e p . This is an 
approximate estimate, since, in the range of crossover at 
temperatures near T q , the local superfluidity and the lo- 
cal long-range order vanish, elementary excitations cease 
to form a nearly ideal gas, and the quantum regime goes 
to the classical regime. 
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